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Abstract 

We report the first observation of the charmless hyperonic B decay, B~^ -^ AAi^"*", using a 140 
fb~^ data sample recorded at the T(4S) resonance with the Belle detector at the KEKB e~^e~ 
coUider. The measured branching fraction is B{B+ -^ AAK+) = (2.9lto:?o ± 0.38) x 10"^. We 
also perform a search for the related decay mode B'^ -^ AAvr'^, but do not find a significant signal. 
We set a 90% confidence-level upper limit of S(i?+ -^ AA7r+) < 2.8 x 10"^. 

PACS numbers: 13.25.Hw, 13.60.Rj 



Charmless hadronic B decays are of great interest since they provide opportunities for 
probing CP violation, as well as for testing our understanding of strong interactions. While 
charmless mesonic modes were first established over ten years ago [1|], B decays to charmless 
baryonic final states such as ppK^ [2], pAvr" [3|, ppK'^, ppn^ and ppK*^ |^ were first seen 
only recently. The branching ratios for these three-body decays are larger than those for two- 
body modes such as B^ -^ pp, for which only upper limits have been reported [5|. Another 
intriguing feature is the threshold peaking behavior commonly observed in the baryon pair 
mass spectrum [2, Id 1^ • Both features were anticipated by theory |6| , but the underlying 



i were anticipated t 
dynamics are still far from understood [T], |8|, |9|, [lO, lllLll2l|. 

In this paper we report the observation of B~^ -^ AAK~^ decay, the first example of 
a charmless B decay to a final state containing two hyperons [l3]. The rate is found to 
be comparable to that of other charmless three-body baryonic modes as well as that for 
B — > (jycjyK 1J| . The invariant mass of the AA system has a prominent near-threshold peak. 



With three strange particles in the final state, the AAK^ mode may complement b — > sss 
dominated mesonic modes such as B ^ (pK'^*^ jlSl Il6l llTI . Ilq , for which the polarization and 
CP asymmetry may be sensitive to new physics. With the three-body final state and self- 
analyzed polarization information from the A decay p, Isl, |lO|, |ll | , the B^ — > AAK^ process 



can be used to probe not only CP violation, but T (time reversal symmetry) violation as 
well. 

We use a data sample of 140 fb^^ integrated luminosity, consisting of 152 million BB 
pairs with no accompanying particles, collected by the Belle detector at the KEKB asym- 
metric energy e+e~ (3.5 on 8 GeV) collider [l^. The Belle detector is a large solid angle 
magnetic spectrometer that consists of a three- layer silicon vertex detector (SVD), a 50- 
layer central drift chamber (CDC), an array of aerogel threshold Cerenkov counters (ACC), 
a barrel-like arrangement of time-of-flight scintillation counters (TOF), and an electromag- 
netic calorimeter comprised of CsI(Tl) crystals located inside a super-conducting solenoid 
coil that provides a 1.5 T magnetic field. An iron flux-return located outside of the coil 
is instrumented to detect K^ mesons and to identify muons. The detector is described in 
detail elsewhere J20|. 

Since KEKB operates with a center-of-mass energy at the T(4S) resonance, which decays 
into a BB pair, one can use the following two kinematic variables to identify the recon- 



structed B meson candidates: the beam constrained mass, Mbc = yE^^^^ ~ Pb-, and the 
energy difference, ^E = Eb — -Ebeam, where -Ebeam, Vb^ and Eb are the beam energy, the 
momentum, and energy of the reconstructed B meson in the T(4S) rest frame, respectively. 
The Mbc resolution of about 3 MeV/c^ is dominated by the beam energy spread. The AE 
resolution for B^ —>■ AAi^"^ ranges from 12 MeV to 17 MeV, depending on M^^. 

The event selection criteria are based on the information obtained from the tracking 
system (SVD+CDC) and the hadron identification system (CDC+ACC+TOF), and are 
optimized using Monte Carlo (MC) simulated event samples. 

All primary charged tracks are required to satisfy track quality criteria based on the 
track impact parameters relative to the interaction point (IP). The deviations from the IP 
position are required to be within ±0.3 cm in the transverse (x-y) plane, and within ±3 cm 
in the z direction, where the z axis is the opposite of the positron beam direction. 

Primary kaon and pion candidates are selected based on K/n likelihood functions ob- 
tained from the hadron identification system. To identify kaons/pions, we require the like- 
lihood ratio Lk(tt)/{Lk + L.„) to be greater than 0.6. For kaons(pions), this requirement 
has an efficiency of 86%(89%) and a pion(kaon) misidentification probability of 8%(10%). 
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FIG. 1: AE and Mbc distributions of B+ -^ AAK+ candidates for M^^ < 2.85 GeV/c^ The solid 
and dashed curves represent the fit results and the signal, respectively; the dotted curve shows the 
background contribution. 



A candidates are reconstructed via the pir decay channel using the method described in 
ref. 0. 

The dominant background for the rare decay modes reported here is from e~^e~' —* qq 
continuum processes (where q = u,d,s,c). The background from generic B decays and 
known baryonic B decays is negligible. This is confirmed using an off-resonance data set 
(10 fb~^) taken 60 MeV below the T(4S) and MC samples of generic B decay, 150 million 
continuum events and known baryonic B decays. In the T(4S) rest frame, continuum events 
tend to be jet-like while BB events tend to be spherical. We follow the scheme defined 
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that combines 7 shape variables to form a Fisher discriminant |22| in order to optimize 
continuum background suppression. The variables used have almost no correlation with Mt,c 
and AE. Probability density functions (PDFs) for the Fisher discriminant and the cosine 
of the angle between the B flight direction and the beam direction in the T(4S) frame are 
combined to form the signal (background) likelihood £s(b)- We require the likelihood ratio 
TZ = Cs/i'C.s + Ch) to be greater than 0.4; this suppresses about 73% of the background while 
retaining 88% of the signal. The selection point is determined by optimizing S/^/S + B, 
where s and b denote the number of signal and background; here a signal branching fraction 
of 4 X 10~^ is assumed. We also require only one candidate per event. In the case of multiple 
B candidates (about 2.6% of the events), we choose the candidate with the highest TZ value. 
The signal PDFs are determined using the signal MC simulation; the background PDFs are 
obtained from the data sideband events with 5.2 GeV/c^ < Mbc < 5.26 GeV/c^ and 0.1 
GeV< \AE\ < 0.3 GeV. 

To ensure the decay is charmless, we exclude 2.85 GeV/c^ < M^^^ < 3.128 GeV/c^ and 
3.315 GeV/c^ < Af^^ < 3.735 GeV/c^ regions to remove contributions from J/^ip, rjc, ip' and 
Xco,i,2 in order to extract the three body decay branching fraction. 

We perform an unbinned extended maximum likelihood fit to the events with —0.15 
GeV< AE < 0.3 GeV and Mbc > 5.2 GeV/c^ to estimate signal yields. The extended 
likelihood function C is 
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l[[NsPsiM^c,,AE,) + NtPtiM^c„,AE,)l 



i=l 



where Ps{Pb) is the signal (background) PDF and Ns{Ni)) denotes the number of sig- 
nal (background) candidates. The signal PDF is the product of a Gaussian function, which 
represents Mbc, and a double Gaussian for AE. The means and the widths of the signal 



PDFs are determined by MC simulation. Differences between data and MC are corrected 
by the B~ — >■ D^n~ and D° -^ K'tt^tt'tt^ control sample. 

We use the parametrization first suggested by the ARGUS collaboration |2^], /(Mbc) oc 

Mbc\Jl - (Mbc/-Ebeam)^exp[-^(1 - (Mbc/-Ebeam)^)] , to modcl the background Mbc distri- 
bution and a 2nd order polynomial for the background AE shape. We perform a two- 
dimensional unbinned fit to the AE-M\yc distribution, floating the signal and background 
normalizations as well as the background shape parameters. 

The Mbc distribution (with \AE\ < 0.05 GeV) and the AE distribution (with Mbc > 5.27 
GeV/c^) for the region Mj^^ < 2.85 GeV/c^ {i.e. below charmonium threshold) are shown 
in Fig. ^ with fit results overlaid. The two-dimensional unbinned fit gives a signal yield of 
22.9l4;8 with a statistical significance of 7.4 standard deviations. The significance is defined 

, where Lq and L^ 



as J—2\n{Lo/Ljnax), where Lq and Lmax are the likelihood values returned by the fits with 
signal yield fixed at zero and floating, respectively 24 . 

We fit the signal yield in bins of M^^ and the result as a function of AA mass is 
shown in Fig^2L_The observed mass distribution peaks at low AA mass, similar to those 
observed in [^, y, yj. Since the decay is not uniform in phase space, we calculate the 
partial branching fraction for each M^;y bin with the corresponding detection efficiency 
determined from a large phase space MC sample and an additional special MC sample 
with a M^^ peak near threshold |,25|. The AA invariant mass spectrum for the events in 
the B+ -^ AAK+ signal region {\AE\ < 0.05 GeV and Mbc > 5.27 GeV/c^) with 2.85 
GeV/c^ < M/^x < 3.15 GeV/c^ is shown in the inset of Fig. |21 A clear J/ijj signal is 
evident. The results of the fits along with the efficiencies and the partial branching frac- 
tions are given in Table EJ We sum the partial branching fractions in Table Q] to obtain 
B{B+ -^ AAK+) = (2.91 to:™ (stat) ± 0.38 (syst)) x 10-^ with a statistical significance of 
5.1 standard deviations. 

The systematic uncertainty in particle selection is studied using high statistics control 
samples. Kaon/pion identification is studied with a D*~^ —>■ D^tt~^, D^ —>■ K^tt^ sample. The 
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FIG. 2: Fitted yield divided by bin size for B^ 



AAK^ as a function of M, 



AA- 



The shaded 



distribution is from a phase space MC simulation with the area normalized to the signal yield. 
Note that the charmonium veto has been applied. The inset shows the AA mass spectrum for the 
ijc and J/V' signal regions. 



TABLE I: Results of the AE — My,c fit, detection efficiencies (e), and branching fractions (B) 
with statistical errors in bins of -/Vf^j^ after the charmonium veto has been applied. The fit allows 
the yields to fiuctuate negative. Note that the yields are consistent with zero above charmonium 
threshold. 
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tracking efficiency is studied with a D* sample, using both full and partial reconstruction. 
Based on these studies, we assign a 7.8% error for the tracking efficiency and 0.6% for 
kaon/pion identification. 

For A reconstruction we have an additional error on the efficiency for off-IP tracks recon- 
struction, determined from the difference of A proper time distributions for data and MC 
simulation. For the four tracks from the AA pair this error amounts to 6.1%. By studying 
the A -^ pn^ sample we assign an error of 1% for each identified proton. There is also a 1% 
error for each A mass selection and a 0.7% error for each A vertex selection J5|. Summing 
the correlated errors for A and A reconstruction, we obtain a systematic error of 6.9% for 
both A's. 

Continuum suppression is studied with the topologically similar B^ -^ D^n^, D° — >■ 
K~7T~^7T~^7T^ sample. By changing the selection criteria on TZ in the interval - 0.4, the 
efficiencies of data and MC differ by 3%. 

The systematic uncertainty from fitting is 4.9%, which is studied by varying the pa- 
rameters of the signal and background PDFs by ±lcr. The MC statistical uncertainty and 
modelling with six M^^ bins contributes a 5.0% error(obtained by changing the M^^ bin 
size). The error on the number of total BB pairs is determined to be 0.7%. The error from 
the sub-decay branching fraction of A — >■ p7i~ is 0.8% |2J]. 

We sum the correlated errors linearly and then combine the result with the uncorrelated 
ones in quadrature. The total systematic error is 13.0%. 

We perform a cross-check of this analysis by using the charmonium- veto events. We 
measure B{B^ — >■ J/ipK^) by following the same analysis procedure with 3.06 GeV/c^ < 
^AA < 3.14 GeV/c^. The signal yield is 11.4l|;2 candidates with a statistical significance 
of 7.3 standard deviations. The obtained product branching fraction is B{B^ — > J/iIjK~^) x 
B{J/ij -^ AA) = {1.55toil) X 10-^ By using B_{B+ -^ J/i)K+) = (1.01 ±0.05) x lO-^Q, 
our measured branching fraction is B{J /ii -^ AA) = (1.54 to;43± 0.20 ±0.08) x 10"^, which 
is in agreement with the world average value [2^. The third error comes from the uncertainty 
oiB{B+ ^ J/ijK+). 

We also search for the decay mode B^ -^ AAvr^, which is an example of a 6 — > uud 
process with ss popping. The background from B^ —>■ AAK^ is negligible. We perform a 
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FIG. 3: A^ and Mk^ distributions for B' 



AAvr^ candidates. The solid curve is the fit result. 



two-dimensional unbinned likelihood fit to the AE'-M^c distribution using the same analysis 
procedure. No significant signal is found. The Mbc and AE distributions with fit projections 
are shown in Fig|21 We use the fit results to estimate the expected background, and compare 
this with the observed number of events in the signal region in order to set an upper limit 
on the yield at the 90% confidence level J26|. The estimated background is 37.5 ± 1.0, the 
number of observed events is 41, the systematic uncertainty is 15%, and the upper limit 
yield is 21.7. The efficiency, estimated from the phase space MC, is found to be 5.05%. The 
90% confidence level upper limit for the branching fraction is B{B^ — >■ AAtt^) < 2.8 x 10~^. 

In summary, we have performed a search for the rare baryonic decays B~^ -^ AAK~^ 
and AAtt"*" with 152 million BB events. A clear signal is seen in the B~^ -^ AAK^ mode, 
where we measure a branching fraction of B{B~^ -^ AAK~^) = (2.9llo:7o =t 0.38) x 10~^, 
which is comparable to B^ -^ ppK~^ and B^ -^ pAn^ . The observed Mj^j;^ spectrum 
peaks toward the threshold as in the above mentioned modes. This measurement is the 
first observation of a i? meson decay to a hyperon pair through a 6 — > sss process. The 
B~^ — > AA7r+ mode is not statistically significant, and we set the 90% confidence level upper 
limit B{B+ -^ AA7r+) < 2.8 x 10"*^. 
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